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Abstract. Orthogonal polynomials of degree n with respect to the weight 
function Wfx(x) = (1 — ||x|| 2 )' 1 on the unit ball in M. d are known to satisfy the 
partial differential equation 

- (x, V) 2 - (2fi + d)(x, V>] P = -n{n + 2/i + d)P 

for fi > — 1. The singular case of fi = —1, —2, ... is studied in this paper. Ex- 
plicit polynomial solutions are constructed and the equation for v = — 2, — 3, . . . 
is shown to have complete polynomial solutions if the dimension d is odd. The 
orthogonality of the solution is also discussed. 



1. Introduction 

Let Wfj,(x) = (1 — H^H 2 )' 1 . For fi > — 1, is integrable on the unit ball 
B d := {x : \\x\\ < 1} and orthogonal polynomials with respect to W„ are well 
defined. We denote by V^iyV^) the space of orthogonal polynomials of degree n 
with respect to W M on B d . It is well known (cf. [5]) that P G V^(W^) satisfies the 
partial differential equation 



(1.1) L^P:= A - (x, V) 2 - (2/i + d)(x, V) 



P = Xi^P, with 



where A denotes the usual Laplacian, V = (d\, . . . ,dd) with di — so that 

(x, V) = Yli=i x i~Sr- ^ n °tber words, the elements of V ? f (W 7 ^) are eigenf unctions of 

Lfj, with eigenvalues A„ . Since U„>oV^(W M ) contains a basis for all polynomials, 
we say that has complete polynomial solutions for fi > — 1. 

The purpose of this paper is to study equation (|1.1|) for \i = —1,-2,.... Our 
study starts from the realization that orthogonal polynomials with respect to the 
inner product 



(/, <?>_i = A / Vf(x) ■ Vg(x)dx + / f{x)g{x)<kj, A > 0, 

JB d JS*- 1 

which were studied recently in [T^], satisfy for fj, = — 1 (see Section 2). To verify 
it, explicit orthogonal structure from [12) is used. It turns out that it is possible 
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to construct explicit polynomial solutions for the equation when —fi G N. 

However, somewhat surprisingly, has a complete set of polynomial solutions 

only when the dimension d is odd, and our construction shows why this is so. A 
natural question is whether the polynomial solutions of (jl.ip remain orthogonal 
for fi = —2, —3, . . . when d is odd. We do not know the answer to this question. 
We are, however, able to show that a family of polynomials closely related to the 
solutions of the equation (jl.ljl for k = — 2 are orthogonal with respect to 

(f,g)_ 2 = \[ Af(x)Ag(x)dx+ f f{x)g{x)dLO, A > 0. 

JB d JS- 1 - 1 

Furthermore, for each n G N, the polynomial solutions of fll.lj) for k — —2 are 
orthogonal with respect to an inner product that is a modification of (•, •)_„ but 
depends on n. 

Historically, only the case d = 1 of such a problem has been thoroughly studied. 
In the case of d = 1, the equation becomes the ordinary differential equation 
satisfied by the Gegenbauer polynomials. In this case, the problem of negative in- 
dices has been studied by several authors, we refer to [TJ[S1 [31 [5] and the references 
therein. For d = 2, equation (jl . ip is classical and can be traced back to Hermite. It 
is one of five second order partial differential equations with polynomial coefficients 
that have orthogonal polynomials solutions with respect to a positive definite mo- 
ment functional, as classified by Krall and Schcffer in [f\. We note that our result 
for [i = — 1 does not conflict with the result of [7], since (/, <j)_ x does not come 
from a moment functional. For d = 2, the case —fx ^ N is considered in [9], where 
the fact that (jl.ip does not have complete polynomial solutions for —fj, G N is also 
observed. Finally, we should mention [8] and [11] for relevant results. 

The paper is organized as follows. We collect background materials in the fol- 
lowing section. The solutions of (jl.ip for fj, G N are constructed in Section 3. 
Orthogonal polynomials with respect to (•, -)_ 2 are discussed in Section 4. 

2. Background and Preliminary 

Let dui denote the Lebesgue measure on := {x : \\x\\ — 1} and denote the 

area of S"* -1 by u>d, '■= fgd-i ^ = 27r d / 2 /r(d/2). We start with the spherical 
harmonics, which are homogeneous polynomials satisfying the equation AP = 0. 
Let Hn denote the space of homogeneous harmonic polynomials of degree n. It is 
well known that 

d (n + d-l\ fn + d-3\ 

The restriction of Y G W„ on S d ^ 1 are the spherical harmonics, which are orthog- 
onal on iS* 1 ^ 1 . Throughout this paper, we use the notation {Y™ : 1 < v < a„} to 
denote an orthonormal basis for 7i d , that is, 

(2.1) — / Y r ;(x')Y™ (x')dw(x') = S MlV 6 n , m , \<^v<o n . 

We will also denote by the space of polynomials of degree n in d variables and 
by V„ the space of homogeneous polynomials of degree n. It is well known that 
dimP^r^- 1 ). 
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For orthogonal polynomials on the unit ball, we start with the inner product 
{f, 9)^ = ^1 f(x)g(x)(l-\\xfrdx, M >-1- 

JB d 

where c M is the normalization constant of W^. Let V„(W)x) denote the space of 
orthogonal polynomials of degree n. A mutually orthogonal basis for Vf^W^) is 
given by ([5]) 

(2.2) P£ v <Ws,x) = P^ n - 23+ ^\2\\x\\ 2 - l)Y?-*(x), < j < n/2, 

where P^ a '^ denotes the Jacobi polynomial of degree j, which is orthogonal with 
respect to (1 — x) a (l + x)@ on [—1,1], and {Y™~ 2 i : 1 < v < <r n -2j} is an orthonor- 
mal basis for Hf l _ 2 j- 

We will need orthogonal polynomials for two other inner products on the ball. 
One is (•, in the introduction, which we normalize as 

(/, g)_, := — f V/(x) • Vg{x)dx + — f f(x)g(x)du d , A > 0, 

^>d J B d u d JS 4 - 1 

so that (1, 1) * = 1. Let V^(W_i) denote the space of orthogonal polynomials with 
respect to (1, Then, as shown in |12j . 

(2.3) V^(W_i) = H d n e (i - ||s|| 2 )v*_ 2 (wi). 

The second inner product is motivated by an application in numerical solution of 
Poisson equations (see [4]) and it is defined by 

(2.4) (/, g) A =a d f A[(l - \\x\\ 2 )f(x)} A[(l - \\xf)g{x)]dx, 

JB d 

where a d = 1 / (4d 2 vol(B d )) so that (1, 1) A = 1. Let V^(A) denote the space of 
orthogonal polynomials with respect to (•, •) » . Then, as shown in [llj . 

(2.5) V*(A) = Hi © (1 - \\x\\ 2 )V*_2(W 2 ), 

We refer to [TTJ [T2] for further properties of orthogonal polynomials with respect 
to these two inner products. Note that an explicit orthonormal basis for either 
orthogonal polynomial subspace can be easily deduced from (|2.2[) . 

In the above and throughout this paper, we write g(x)V n = {g(x)p(x) : p 6 
V n }, where V n is either the subspace H^ of harmonic polynomials or the subspace 
Vn(Wfj,) of orthogonal polynomials. Furthermore, if n < 0, then we define H„ = 
and Vi = 0. 

3. Polynomial Solutions of the Partial Differential Equation 
We start with two lemmas on the operator L^f. 

Lemma 3.1. For /x € K. 

(3.1) L„[(l - \\x\\ 2 )P] = 4(/i + 1)P + (1 - \\x\\ 2 ) [-(4(/x + 1) + 2d)P + L M+2 P] . 
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Proof. The proof follows from a straightforward computation. The essential ingre- 
dients are 

A[(l - \\x\\ 2 )P] =-2dP- A{x, V)P + (1 - \\x\ 2 )P, 
(x, V)[(l - \\x\\ 2 )P] = - 2P + 2(1 - |M| 2 )P + (1 - \\x\\ 2 )(x, V)P + (1 - \\x\ 2 )P, 
(x, V} 2 [(1 - |M| 2 )P] = - 4P - 4(x, V)P + 4(1 - |M| 2 )P 

+ 4(lHM| 2 )(z,V)P+(l-|M 2 )<z,V) 2 P, 
as can be easily verified. □ 

Lemma 3.2. For all fie R and n>0, L^Y = A^Y for Y e H d ri . 

Proof. If Y £ U d n , then AY = 0. Since Y is homogeneous, Euler's formula shows 
that (x, V)Y — nY. The stated result follows immediately from these two facts. □ 

Our first result shows that orthogonal polynomials with respect to (•, satisfy 
the equation (jTTTJ) with \i = —1. 

Theorem 3.3. Elements of V d {W-\) satisfy L-±P = l 'P. In particular, the 
eigenf unctions of the operator L_i consist of a complete polynomial basis. 



Proof. Setting fj, = — 1, equation (|3.1|) shows that 

L_i [(l-\\x\\ 2 )P}=(l-\\x\\ 2 )[-2dP + L 1 P}. 

Recall that elements of V„_ 2 (Wi) satisfy the equation (|1.1[) with /i = 1. It follows 
that for P e V£_ a (Wi), 

L_i [(1 - |M| 2 )P] -(1 - |M| 2 ) \-2dP - (n - 2)(n + 2// + d - 2)P] 

= -n(n + d-2)(l- |M| 2 )P. 

By (|2.3p . this and Lemma T3.2I show that all elements in V„(W_i) satisfy (jl.ip for 
H = -1. □ 

For <i = 2, a complete characterization of second order linear partial differential 
equations that have polynomial coefficients and orthogonal polynomials as eigen- 
functions is given in [TJ, where the orthogonality is given in terms of a moment 
functional. If the moment functional is positive definite, then there are only five 
families of equations, the equation p.l[) is one of them when fj, > — 1. It should be 
remarked that Theorem 13.31 does not contradict the characterization of [7J , since 
the inner product (•, cannot be expressed as a moment functional. Indeed, it 
is easy to see that (xf,g)_ 1 7^ (f,xg)_ 1 in general. 

The above theorem prompts us to consider solutions of (jl.ip for fi = —2, —3, . . . 
and leads to our second theorem. For ji = —k, k — 2, 3, . . ., we define 

' ni_ 2 A\j{i-\\x\\ 2 ) k v d n _ 2k {w h ), 



5>2„(i-IM 



u*(W- k ) -.= n d n u 

where, for 1 < j < k — 1, 

3.2 c#„ := — — 0<v<j, 

1 ' hV ulij-uy.i-k + lUn-j-k + Dj' ~ ~ Jl 

in which (a) m := a(a + 1) . . . (a + m — 1) is the shifted factorial. We note that 
djj = 1. It is easy to see that a™„ are well defined if n — j — fc + ^ + | is not zero. 
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Theorem 3.4. If fj, = —k and k — 2,3,..., then the polynomials in U^W-k) 
satisfy equation that is, L_^P = Xn k ^ P for P 6 U^(W-k)- Furthermore, 

dimU% = dim V% if n - 2k - 1 + f ^ 0. 

In particular, the operator has a complete polynomial basis of eigenf unctions 
if the dimension d is odd. 

Proof. We observe that the relation (|3.ip can be used recursively, which yields 
(3.3) L,{l~\\xf) k =Ak^ + k){l-\\xf) k - 1 

+ (1 - ||x|| 2 ) fe [-(4*c(/i + k) + 2kd) + Lp+tk] , k e N, 

as can be easily verified by induction. Hence, for Q G V^_ 2fe (Wfe), which satisfies 
LkQ = A„_ 2 fc < 3> we nave 

L_ fc [(l - ||x|| 2 ) fe Q(x)] = (I - ||x|| 2 ) fe \-2kdQ(x) + L k Q(x)] 

= (1 - ||a;|| 2 ) fe \-2kd - (n - 2*)(n - 2fe + 2k + d)} Q(x) 

= -n(n -2k + d)(l - ||a;|| 2 ) fe Q(x), 

so that elements in (I - |M| 2 ) fe V^_ 2fe (W0c) satisfy L- k P = Ai~ fe ' P. Furthermore, 

by Lemma l3~2l elements of are also eigenfunctions of L_fc with eigenvalue A« . 
Let now 

j-i 

Z,0r) = (1 - |]ar|| 2 )iy(ar) + ]T - INI 2 ) 1 ^), y g ^_ 2j . 

!/=0 

We want to show that Zj satisfies the same equation. Applying (|3 - 3[) shows that 

j 

L-kZjix) =-n(n-2k + d)(l - \\xf) j Y [x) + [ 4 (" + !)(- fc + v + l ) a t»+i 

!/=0 

-(n - 2j + 2u){n - 2j + 2v-2k + 2)o!?„] (1 - \\x\\ 2 ) v Y(x), 

so that the equation L-^Zj = \h k ^Zj implies, upon comparing the coefficients, 
that a^l satisfy the recursive relation 

in) («) -2(v+l){-k+v+l) 

a) ' = a) 



3> j ' v+1 (j - i/)(2n - 2j + 2v - 2k + d) ' 

Setting — 1, we can then determine recursively. The result is (|3.2p . This 

shows that all elements in l4£(W—k) satisfy the equation L-kP = A« fe ^P. 
If a™ v are all finite, then it follows from dimP^ = dim 7-^ + dim'P„_ 2 that 

dimW*(W_ fc ) = dimtt* + . . . + dim«£_ 2fc+2 + dimP*_ 2fc = dim 7^. 

The explicit formula (|3.2p shows that a™^ is not finite only if n — j — k + v + d/2 is 
zero, which can happen only if d is even and n < j + k — v — d/2. Since < v < j — 1 
and j < k — 1, this could happen only if n — 2fc — I + d/2 is zero or a negative 
integer. □ 
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Some remarks are in order. If d is even then, as shown in the proof, U^(W-k) is 
not well defined if n — j — k + d/2 = 0. The first time this happens is when d = 2 
and k = 2. As 1 < j < k — 1, this happens when j = 1 and n = 2. In fact, for 
fc = 2, we have 



Hi U 



(1 



1 



^- 2 U(l-||x| 2 ) 2 V2_ 4 (^ 2 ), 



so that U^{W-2) is not well defined. Furthermore, since A 2 = 0, H2 = spanja; 2 — 
x|, 2x1X2} and L^l = 0, we see that the linearly independent eigenfunctions with 
respect to the eigenvalue A 2 2 ' = are l,a; 2 — x\,2x\X2- However, 1 is also the 
eigenfunction of Aq , we see that L_ 2 does not have a complete polynomial basis 
as eigenfunctions and it misses exactly one polynomial of degree 2. 

The same analysis also works for k > 2. For example, for k = 4 and d = 2, 



Ul{W. 



Hi u 



u 



u 



(1 
(1 
(1 

+- 



N 2 



|2\2 



|2\3 



3 


n — 


4 


4 




n — 


4 



i_/2 



n - 4 
G 



(i-IMI 2 ) + 
(i-INI 2 ) 2 h 



ri n-6 



(n - 4)(n - 5) 
6 

" (n-4)(n-5) 

u(i-]N 2 ) 2 v^ 



ri n-4 

(i-IMI 2 ) 

-8(^2), 



(n-4)(n-5)(n-6)_ 
which shows clearly that the operator L_4 has a complete basis of polynomials 
of degree exactly n as eigenfunctions with respect to the eigenvalue A„ only if 
n > 6. 

In general, for even dimension d, the operator L_k has a complete basis of 
polynomials of degree exactly n as eigenfunctions only if n > 2k — 1 + d/2. As 
a consequence, the eigenfunctions of the operator do not contain a complete 
basis of polynomials if d is even. 



4. Orthogonal polynomials associated with Laplacian 

For all fi € K, the function W" M is a symmetric factor for the differential operator 
Lfj as it is easy to verify that L^W^ = (W^L^)* with 

{WpLJ* := AW,, + (d + 2 M )(V,x)W M - (V,x) 2 ^. 

The verification is formal and holds for all The polynomial eigenfunctions of 

I/ M are orthogonal with respect to W M for /i > — 1. According to Theorem 13.31 the 
polynomial eigenfunctions of L_i form a complete orthogonal basis with respect to 
(•, It is natural to ask if the same holds true for L-k with fc = 2, 3, . . . for d 
being odd. Below we consider the case k — 2. 

Recall that in the case of k = 2, the eigenfunctions of L_2 are given by 

n d n _ 2 u(i-\\ x ffvtM. 

We first consider the following closely related class of functions 

(4.2) V*(W-t) := H d n U(l- ||x|| 2 )^_ 2 U (1 - ||x|| 2 ) 2 Vl 4 (^ 2 ). 



(4.1) U*(W- 2 )=Hiu 



(1- w 5 



n-3+l 
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The difference between V d (W-2) and Vf^W-2) are only in the middle term. Notice 
that V„(W-2) are defined for all n and d. It turns out that it is the space of 
orthogonal polynomials with respect to the inner product 

(/, g)_ 2 := — f Af(x)Ag(x)dx + — f f(x)g(x)dw, A > 0, 

JB d ^d J gd-i 

which is normalized so that (1, 1)_ 2 = 1- It is easy to see that this inner product 
is well-defined and it is positive definite. 

Theorem 4.1. The elements 0/ V^(VF_ 2 ) are orthogonal polynomials with respect 
to (■, ■) _ 2 - Moreover, they contain an orthonormal basis; in other words, 

V d n {W- 2 ) = K © (1 - \\xf)H d n _2 © (1 " IM| 2 ) 2 V^ 4 (W/ 2 ). 

Proof. First we make an observation that, by (|2.5[) . 

V„ d (M/_ 2 )=^U(l-||xf)V r t 2 (A). 

We then prove that each of the two components in the right hand side of the 
above expression are orthogonal to lower degree polynomials and they are mutually 
orthogonal with respect to (■, -)_ 2 - We consider the two cases separately. 

Case 1. Let P = Y n g U d n . Since AP = 0, {P,g)_ 2 = ^ J sd -i f(x)g{x)dto. If 
g g II^_ 1S then (P,g)_ 2 = by the orthogonality of Y n . If g g (1 - ||ir|| 2 )V^_ 2 (A), 
then the restriction of g on S^ -1 is zero, so that (P,g)_ 2 = 0. 

Case 2. Let P(x) = (1 - \\x\\ 2 )Q n ^2{x), where Q„_ 2 6 V d _ 2 {A). Then the 
restriction of P on is zero, so that 

(P,g)_ 2 = - [ A[(l - \\x\\ 2 )Q n -2{x)]Ag{x)dx. 

Let g £ V^W-2) and m < n. If g g W m , then Ag = so that {P,g)_ 2 = 0. 
If G (1 - |M| 2 )V* _ 2 (A), then 0(3) = (1 - ||a:|| a )/*(a;) for h g V^_ 2 (A). Recall 
the definition of {f,g) A in (|2.4[) . we then have (P,g)_ 2 = (Q n _ 2 ,/i) A = by the 
orthogonality of Q n -2 6 V^_ 2 (A). □ 

Recall that V^(IF_ 2 ) and l^*(W_a) differ only at the middle component. We do 
not know if there is an inner product with respect to which £/„(W_ 2 ) is orthogonal 
for all n. For each fixed n <E N, however, elements of W„(W_2) are orthogonal to 
lower degree polynomials with respect to the following inner product: 

(f,9)n-=— [ Af(x)Ag(x)dx + — [ f{x)g(x)doj 

Ud J B d w d J 3 d - 1 

+ -[ T [r n - A ~ d m] T [r n - 4 - d 9(x)] du, 
Lu d J S d-i dr dr 

where -£ is the normal derivative, A > and fx > 0. It is easy to see that (/, g)* n is 
positive definite for f,gE H d . 

Theorem 4.2. The elements of V^(VF_ 2 ) are orthogonal polynomials with respect 
to the inner product (f,g) n . 
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Proof. Notice that the first two terms in the (/, g)* n are exactly (/, g)_ 2 - For x G M. d , 
write x = rx', r > and x' G S . We assume that <? € n^_ 1 below. If P = Y n G 
7Y^, then ^ [r"~ 4 - d P(x)] = (2n-4-d)F„(a;') on S"*" 1 . Hence (P, #)* = as AP = 
and P G CH%. If P G (1 - ||z|| 2 ) 2 V^_ 4 (W 2 ), then £ [r"- 4 - d P(a;)] = on S^ 1 
because of the (1 - r 2 ) 2 term, so that {P,g)* n = (Ps)_ 2 = 0. Thus, by (pL"T]) . we 

are left with the case P(:r) = (1 — ||a;|j 2 ) + n _ 3 1 +d / 2 Y n ^ 2 (x), where Y n - 2 G W^_ 2 - 

Since Y n _ 2 is homogeneous, Y n _ 2 {x) = r n Y n _2(x'), a simple computation shows 
that £ [r n ~ 4 - d P(x)] = 0. Consequently we again have (P, g)* n = (P, g) _ 2 =0. □ 

These results indicate that there is a difference between the case k = — 1 and 
k = —2,-3,.... Not only the differential operator L- k has complete polynomial 
eigenfunctions only when d is odd in the latter case, the eigenfunctions also may 
no longer be orthogonal with respect to an inner product for all n. 

Another interesting question is if the main terms of U d (W-k) are orthogonal for 
k > 3. In other words, if the spaces 



V« n {W- k ) := n d n U |J (1 - \\xfyn d n _ 23 U (1 - \\xf) k V d n _ 2k (W k ) 




are orthogonal with respect to a positive definite inner product, which will likely 
involve higher order derivatives, as indicated by our result for k = 1 and k = 2. At 
this point we do not see a way to find such an inner product. 

Acknowledgment. The authors thank Jeong Keun Lee for communicating and dis- 
cussing the results in [9]. 
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